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Abstract 


This report documents the cavitation assessment for the Space Shuttle 
Main Engine (SSME) high pressure oxidizer turbopump (HPOTP). In order to do 
this, a model of the flow through the pump had to be developed. Initially, a 
computational procedure was developed to analyze the flow through the inlet 
casing including the prediction of wakes downstream of the casing vanes. From 
these flow calculations, cavitation patterns on the inducer blades were 
approximated and hence the damage rate estimated. The model correlates the 
heavy damage on the housing and over the inducer with unsteady blade surface 
cavitation. The unsteady blade surface cavitation is due to the large 
incidence changes caused by the wakes of the upstream vanes. Very high 
cavitation damage rates are associated with this type of cavitation. As a 
result of the computational procedure, a set of design recommendations for 
reducing the unsteady cavitation are given. These include removing the set of 
vanes closest to the inducer and modifying the remaining vanes. 
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l ♦ Introduction 

The destructive action caused by cavitation has long been a problem of 
practical importance with fluid machinery. At the turn of the century, severe 
erosion of the propellers on the ships Lusitania and Mauretania [1] prompted 
the creation of a special commission by the British Admiralty to investigate 
the problem. It was determined chat the damage was the result of repeated 
"hydraulic blows" accompanying the collapse of cavities. 

After nearly 100 years of research in this area, the problem of 
predicting damage on a prototype is still unsolved. This is not due to a lack 
of work in this area but rather due to the fact that cavitation damage 
involves both solid mechanics and the fluid mechanics of two-phase flows, and 
thus is inherently difficult. As nottsd in a recent article by Hutton (2] 
there has been some success in empirically predicting when and where 
cavitation may begin, but there are difficulties in predicting how serious its 
damage effects might be. The literature on cavitation damage is very 
extensive and as noted by Hutton [2] "provides absorbing reading for those 
with a morbid interest in engineering failures". 

The problems associated with the occurrence of severe cavitation damage 
in fluid handling machinery are primarily dependent on the degree and type of 
cavitation. The majority of all fluid machinery can operate with some limited 
degree of cavitation and still provide satsifactory performance although 
damage can still occur if operated over longer periods of time. However, 
large amounts of cavitation not only causes severe cavitation damage, but the 
accompanying rapid formation and collapse of cavities alters the local 
pressure field acting on each blade segment. As a result, net axial and 
radial thrusts become random, reaching high magnitudes at unpredictable 
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frequencies. The unsteady forces generated in this manner cause severe wear 
and damage to bearings and seals. 

The practical problem addressed in this report is the reduction of 
cavitation damage in the Space Shuttle Main Engine (SSMf high pressure 
oxidizer turbopump (HPOTP) and the improvement of its life expectancy. 
Experience with cavitation damage suggests that this Is a difficult problem 
because it involves (1) a prediction of the pressure field in the turbopump, 

(2) a prediction of the cavitation type and extent, (3) the determination of 
cavitation impact pressures, and (4) an analvsis of the interaction between 
cavitation Impact pressure and the material response. All four of these 
issues nust be addressed in order to predict a cavitation damage rate and 
hence the turbopump life expectancy. This problem is further complicated when 
one realizes that the turhopump will operate off-design for extended periods 
of time. Thus, before developing a damage model, some assumptions nust be 
made from a detail assessment of the existing turbopurap data. This assessment 
includes hardware geometry, hardware evidence of cavitation damage, operating 
conditions and life requirements. 

The characteristics of the SSME are given in Table 1 and the SSME 
powerhead component arrangement showing the HPOTP is illustrated in Figure 1. 
The HPOTP is shown in Figure 2 and its operating conditions are listed in 
Table 2. 

Examination of the turbopump hardware shows that cavitation damage 
occurred at the following locations: 

1. Housing (casing) upstream and alongside the inducer, 


2 . 


Inducer blades near the tip, 
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3* Impeller near the partial blades, and 

4. Impeller blades at the exit diameter. 

Damage on the impeller is relatively light; however, the damage associated 
with the Inducer is very severe. 

figure 3 shows schematically that the damage on the inducer occurs in an 
area located at the blade tip. The most severe damage, however, occurs on the 
housing, downstream of the stator vanes and near the leading edge of the 
inducer. This damage appears at discrete spots as can be noted in Figure 4 
for the left-hand pump housing and in Figure 5 for the right-hand pump 
housing. In addition, a band of cavitation damage occurs on the housing near 
the trailing edge of the inducer as can be noted in Figure 6. 

The deep-pitting cavitation damage occurring on the housing significantly 
reduces the life expectancy of the turbopump. A depth of cavitation damage of 
0.05 inches in only 30 minutes of operation has been documented by Rocketdyne 
and Figure 7 shows the characteristics of a typical pit. This time is 
significantly less than the goal of 7 .5 hours for the life expectancy of the 
turbopumo (see Table 1). 

High damage rates can be expected in this area of the turbopump because 
of high inducer tip speed which is in excess of 600 ft/sec. Many 
investigators such as ’(napp [31 and Stinebring [4] have shown that the damage 
rate is proportional to the sixth power of the relative velocity. 

The existence cl deep pitting suggests that bubble collapse is occurring 
predominantly at these locations. This is possible if the amount of 
cavitation occurring on the inducer varies circumferentially, and could result 
from a changing flow incidence angle. 

In this report, a model is developed which correlates the heavy damage on 
the housing and over the inducer blades with unsteady blade surface 
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cavitation. Initially, a flow calculation procedure waa developed that gives 
the inlet flow distribution at the Inducer Including the wakes from the 
upstream vanes. The cavitation type and extent were predicted using the inlet 
data, flow conditions and blade geometry. The cavity volume shed was then 
used to predict the cavitation damage rate. Finally , recommendations for a 
redesign are given that should significantly reduce the cavitation damage 
rate, 

2. Analysis of Existing Design 

2.1 Inviscid Flow Calculations 

The flow through the High Pressure Oxidizer Turbopump (HPOTP) is very 
complicated and the design of its components is critical. It consists of five 
blade rows, three of which are stationary, and the flow is turned four times 
in a very short distance. This compact geometry is illustrated in Figure 2. 
The blade rows are identified as (l) the inlet casing ''anes, (2) the stator 
vanes, (3) the inducer, (4) the centrifugal impeller, and (5) the exit 
diffuser vanes. The flow analysis begins with an examination of the first two 
blade rows. 

The inlet casing is a very difficult component to analyze. Flow enters 
through a pipe and is split in half to form a double-suction entrance. The 
flow is then divided again in each side inlet (left and right), and is thus 
symmetric about a vertical axis, but each side inlet is asymmetric about the 
axis of rotation. Only a three-dimensional (3-D) code would be able to give a 
true prediction of the flow in this element. While there may exist 3-D codes 
to solve for the inlet casing alone, it is doubtful that a code exists at this 
time that would be able to include the stationary vane rows as well. 

Since flow through the inlet casing is predominantly radial for the 
stationary vane rows, a two-dimensional (?.-D) approximation for the inlet 
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casing was chosen. A side view of this geometry is illustrated in Figure 8 
and allows computation of the flow components in the radial (R) and 
circumferential (6) directions only. The flow will be considered to be 
uniform in the axial (Z) direction. 

Based on the design drawings of the HPOTP, a 2-D configuration can be 
constructed for the inlet casing boundary, the inlet casing vanes, and the 
stator vanes which is illustrated in Figure 9. The fifth stator vane at the 
bottom of the inlet has been blended into the casing boundary and provides for 
vertical symmetry of the flow in the inlet. It is still impractical, however, 
to solve for the flow in this configuration. A transformation must be 
performed to convert this radial flow geometry into a linear cascade 
geometry. 

It is possible to solve for the flow in this geometry if all the 
stationary vanes are removed. The resulting singly-connected domain can be 
solved by using a finite element method. A grid-work of triangular elements 
is illustrated in Figure 10. A total of 647 nodes and 1180 elements were used 
in this configuration. 

If the flow is assumed to be incompressible (Mach number below 0.3), the 
continuity equation in cylindrical or polar coordinates can be written as 
follows: 



1^ + V_r 

r 30 r 


0 


This equation can be made an identify by defining a streamf unction: 


( 1 ) 


1 3ij> 
r 30 



( 2 ) 


Thus, one has 
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For irrotational flow (w-0), the Laplace equation results from the definition 
of angular momentum and can be written as 


-wr 


38 


3r 


- V 9 


ar 2 9r " r a0 ; 


(4) 


This equation can be readily solved by even the most basic of finite element 
me thods . 

The boundary conditions for the strearaf unction are a value of 0.0 on the 
right vertical boundary and a value of 1.0 on the outer casing boundary. The 
value of increases linearly from right to left for the top boundary and from 
too to bottom counterclockwise for the inner boundary. This results in a 
uniform normal velocity at the top inlet and at the inner circular exit. This 
does n^t mean that the tangential velocity at these locations will necessarily 
be equal to zero. Because of the geometry of the top inlet, the flow is 
everywhere parallel to the side boundaries and there is nothing capable of 
sustaining a pressure gradient across the channel which would allow a 
tangential component of the flow to exist. As the radius of the inner 
circular boundary is reduced to zero the flow here nust also be uniform, and 
thus no tangential flow component will exist at this boundary. There will 
however exist a substantial Vg component at the radii of interest, that is the 
exit of the inlet casing vanes and the exit of the stator vanes. 

The finite element solution is illustrated in Figure 11. The streamlines 
that are plotted (constant values of ty) represent equal increments of flow. 
Lines that are close together indicate high velocities and lines that are 
further apart indicate lower velocities. It can be seen that the flow is 
indeed uniform at the inlet and the exit. At the indicated radii a 
substantial Vg component of the flow is shown to exist. It will be the task 
of the stationary vanes to remove this tangential velocity. 


-15- 


The solution technique employed to solve for the 2-D configuration with 
vanes is the Modified Douglas Neumann Cascade Program [6] . In this approach 
each body or vane is divided into many small segments and at each is located a 
source/sink of unknown strength. The Neumann boundary condition is used by 
this method and requires that there be no flow normal to any boundary. For 
lifting bodies (those that are not symmetric normal to the flow direction) , 
the Kutta condition is also employed to insure that the flows from both 
surfaces exit the body in parallel at the trailing edge. For a specified flow 
angle (for example the inlet flow angle) a solution is obtained for all of the 
source/sink strengths, the accuracy of which reaches one hunored percent as 
the number of points or segments on the bodies becomes infinite. With these 
values known the velocity can be calculated anywhere in the field. At 
locations far upstream and downstream the calculated velocity will necessarily 
be parallel and uniform. 

The original Douglas Neumann Cascade Program was developed by the Douglas 
Aircraft Corporation [5] and extended the existing solution technique being 
employed for isolated airfoils to cascades of lifting or non-lifting bodies. 
The program has been modified to include cascades with radial flow, either 
inward or outward, and cascades with bodies of varying geometry. Radial flow 
cascades are solved by transforming them into linear cascades using a 
conformal transformation: 

x = in r 

(5) 

y * 9 

Velocities calculated for the linear cascade are converted to the radial-flow 
plane by simply dividing by the radius as follows: 

= y , Vq = “ Vy 


( 6 ) 
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Special considerations must be taken to account for the inlet casing boundary 
since one would not normally consider this element to be a body in a cascade. 

If the inlet casing boundary shown in Figure 9 were transformed using 
Equation (5) the large body with the dotted lines shown in Figure 12 would be 
produced. With the flow proceeding from left to right it can be seen that a 
large portion of the flow (aproximately eleven-twelths) would be blocked by 
the inlet casing, body number 1. This is due to the flow entering the casing 
through the vertical inlet which only accounts for about 30 degrees of the 
circumference at a radius of 12 inches. It is therefore logical and accurate 
to represent the inlet casing boundary as the streamlined body shown in Figure 
12 to facilitate the flow calculation procedure. 

It is important to model the flow in the vertical inlet correctly. In 
Figure 2 the flow is uniform and directed downward toward the axis of 
rotation. At large radii this flow can be approximated as radial. This will 
result in a horizontal portion of the boundary in the transformed plane. 

Figure 12. Uniform linear flow in this portion of the transformed geometry 
most accurately models the vertical section of the inlet casing. This does 
result in boundaries in the 2-D geometry of Figure 9 which are radial at large 
radii. This type of adjustment is necessary for closure of body number 1 in 
the transformed cascade of Figure 12. 

The Modified Douglas Neumann Cascade Program was first run for only one 
body, that of the inlet casing boundary. For the transformed geometry, 159 
points were used to describe this body. The resulting flow solution is 
illustrated in Figure 13. If this figure is overlayed with the finite element 
solution shown in Figure 11, the two can be seen to correlate exactly. Thus 


one method validates the other 
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The thirteen inlet casing vanes and four stator vanes were then added to 
solve for the flow in the 2-D configuration of Figure 9. Although the vanes 
are symmetric about a vertical axis they are each unique when plotted in the 
transformed plane. The inlet flow angle was set equal to zero to reflect this 
symmetry. The number of points used to describe these vanes ranged from 31 
for the stators to over 50 for the largest Inlet casing vane. In all, 780 
points were used to describe this cascade consisting of 18 bodies. The 
spacing used for each body is 2tt such that the entire geometry illustrated in 
Figure 12 is repeated at the top and bottom to form an infinite cascade. It 
took over an hour of CPU time to solve this problem on the VAX 11/782 
computer. 

The flow solution is first illustrated in the transformed plane in Figure 
14. The velocity vectors are plotted based on information at over 300 
off-body points. In Figure 15 the flow solution has been transformed back to 
the radial 2-D configuration. When compared with the solution for the inlet 
casing boundary alone in Figure 13 it can be seen that the vanes have indeed 
taken out a large portion of the turning. However, a substantial amount of 
tangential velocity still exists and this will have a marked effect on the 
inlet flow at the inducer. 

The flow through the entire cross-section of the HPOTP can be analyzed 
using a streamline curvature method (SCM) [7] . In the SCM program used at 
ARL/PSU the effect of each blade row is accounted for by specifying the blade 
advance ratio (flow coefficient) and the particular distribution of V g at the 
trailing edge station. The Vg profile determines the change in angular 
momentum effected by the blade row and the advance ratio computes the change 
in energy if the blade row is rotating. The SCM program can only solve for 
axisymmetric flow, which is essentially 2-D in that no derivatives are allowed 


%. 
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in the circumferential (0) direction. It solves the continuity, momentum and 
energy equations for the real variables: velocities and pressure. In the 

process, the streamline locations are found as well. 

The streamline curvature solution is illustrated in Figure 16, The 
location of the blade rows are highlighted, including the split impeller 
blades in the rotor. The exit diffuser vanes were omitted for this analysis, 
but could have easily been included. The inlet portion, although represented 
correctly, is not axisymraetric and merely serves to feed the flow into the 
axisymmetric stator, inducer and impeller portion of the flow channel. A 
significant circumferential variation in the flow field was previously found, 
since the channel geometry is not axisymmetric until just prior to entering 
the stator, station 14. It can be seen however that a plane passed normal to 
Figure 16 would closely correspond to the two-dimensional configuration used 
in Figure 8. 

The flow solution illustrated in Figure 16 is for the case of no turning 
added by the blade rows, which for the HPOTP means that the stationary vanes 
have taken out the turning imparted by the asymmetric inlet. There is no 
appreciable difference between the location of the streamlines when a uniform 
(2-D) tangential velocity is added across the span of the inlet casing vanes 
and/or the stator. Essentially the SCM solves for the through-flow and this 
must be combined with the Douglas Neumann circumferential solution to obtain 
the 3-D picture of what is happening at the inducer inlet. 

It now becomes necessary to describe the flow along the various 
streamlines in Figure 16 as the flow turns the corner. Five streamlines from 
the SCM solution were selected which divide the flow into four equal portions. 
There are three locations of primary importance in this part of the analysis: 
(1) downstream of the exit of the inlet casing vanes (R - 3.5 in), (2) the 
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exit of the stator vanes (2.13 in R 2.46 in), and (3) the inducer 
measuring station (1.325 in <. R <. 2.37 in). The corresonding stations for the 
SCM solution are shown in Figure 16 as numbers 13, 18 and 22. The flow will 
be computed at these locations as it turnes the corner. 

The first location at a radius of 3.5 inches will illustrate the flow as 
it exits the inlet casing vane row and enters the stator vane row. Here the 
flow is assumed radial and is calculated from the Douglas Neumann solution of 
the vane geometry in Figure 9. Data points are distributed about the 
circumference at every ten degrees. These are spline fit and plotted in 
Figure 17. The velocities are dimensionalized for the design condition at 
100% power level, but could easily be changed to another condition by 
multiplying the ratio of the volume flow rates at these two conditions. The 
zero degree location corresponds to the top of the circle at this radius, just 
downstream of inlet casing vane number 2. Although this is a potential flow 
solution the obvious location of the stator vanes illustrates the blockage 
that results with these two vanes in the flow on this side of the pump. It is 
at this radius that the wakes produced by the inle*. casing vanes will be added 
to the potential or inviscid flow (see Section 2.2). 

In a similar manner the second location illustrates the flow exiting the 
stator vane row. Again the flow is assumed radial and is calculated from the 
Douglas Neumann solution. Although the stator was calculated as a 2-D blade 
shape the actual length of the vane decreases from hub-to-tip. Five radii are 
therefore employed at station 18 to account for this variation. Since the 
stator vane exit radius was fixed in Figure 9, the blockage effect at the tip 
will be more pronounced than that at the hub. This is illustrated in Figures 
18 through 22 for each of the five selected streamlines. 
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For both of these locations the radial velocity is integrated to check on 
the massflow and diraensionalize the velocity components. The average or 
normalized velocity is indicated for each of the radial positions plotted. 

When the integration is performed the stream function is calculated as a 
function of circumferential location. From Equation (2) with r constant the 
streamf unction becomes 

7T 

=* r T V r A8 . (7) 

o 

These values are normalized by the design massflow to produce values of the 
streamf unction from zero to one as theta increases from 0 to tt radians. 

Since constant values of represent streamlines, the location of a 
selected fluid particle can be calculated at each of the two radii. This 
links up the locations of importance as the flow turns the corner. 

The third location at the inducer station where experimental data are 
available is analyzed in a similar fashion with two major exceptions. First, 
the streamline changes from the radial direction to the axial direction which 
results in a different average through-flow velocity with which to normalize 
the flow due to area changes. Second, the law of constant angular momentum 
which is employed in the SCM cannot be used for this portion of the analysis. 
Recall that the SCM solution is valid for axisymmetric flow only. As the 
fluid particle changes radius from that at the stator vane exit to that of the 
inducer inlet, the casing will continue to act on it and produce a variation 
in the angular momentum of the fluid around the circumference. Thus the 
tangential velocity will not always increase as the radius decreases. 

The potential flow field at the inducer station is depicted in Figure 23. 
There is a noticeable indication of blockage at the two stator locations at 
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rad tl towards the tip streamline. This results from the close proximity of 
the stator vanes and the inducer leading edge at the outer tip boundary. It 
is important to note that the massflow falls off towards the bottom of the 
inlet. This would be expected for an inlet with flow entering at the top. 

2.2 Addition of the Viscous Effects to the Invlscid Flow Analysis 

The viscous boundary layer growth on the stationary blade rows will lead 
to flow deviation and wake velocity deficits downstream of the vanes. These 
effects can readily be included in the flow analysis. In the real fluid flow 
case turning vanes will always achieve less turning than that indicated by 
their camber lines. The fluid will not exit at the blade angle, but at a 
lesser angle. In the inviscid analysis the flow is constrained by the Kutta 
condition to exit at the specific trailing edge angle. In order to achieve 
the proper flow condition downstream of the vane rows, the fluid deviation 
must be added to the inviscid calculation. This was accomplished by fixing 
the vane shapes at the leading edge and increasing the stagger angle so as to 
force the flow to exit at the deviated flow angle. This is illustrated in 
Figure 24. Thus the flow downstream of the vanes will be at the corrected 
angle. 

The deviation angle was estimated from the blade geometry and the blade 
incidence angle obtained from the inviscid flow analysis. The deviation can 
be expressed as 

5 - 5 0 + A6* - AS* (7) 

where 6 0 is the deviation angle due to camber, A 6* is the deviation angle due 


to vane thickness and A6 is the deviation angle due to flow acceleration. 
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The deviations due to thickness and camber were obtained by utilizing 
empirical data which have been collected by Lieblein [8] . The work by 
Lakshmlnarayana [9] was used to estimate the deviation due to flow 
acceleration. 

It is important to note that cascade relationships were used to estimate 
the deviation angles for the vanes inside the casing. However since all of 
the vanes in the HPOTP are of different geometry, the calculated deviation 
angles are only approximate. 

The resulting Douglas Neumann solution shown in Figure 25 indicates an 
even greater disparity of the mass flow from top to bottom. Also due to 
viscous effects the casing vanes will fail to reduce the tangential velocity 
even as much as shown in Figure 15. 

The calculated wakes are also added to the flow at the exit to the stator 
vanes. Effectively the potential velocity distribution is multiplied by a 
wake deficit factor which will be greater than one between vanes and much less 
than one at the vane trailing edge location. The wakes will have different 
widths and depths due to the process of accounting for the wake diffusion 
prior to reaching the inducer leading edge. From Figure 16 it is obvious that 
the wakes at the hub streamline will be able to diffuse (or decay) 
substantially prior to entering the inducer, whereas the wakes at the tip 
streamline will remain narrow and deep. The wakes are diffused by introducing 
the streamwise distance between the vane exit station and the inducer station. 
This wake deficit factor was obtained in the following manner. 

It was assumed that the leading edge of the inducer is located 
sufficiently downstream of the casing and stator vanes that the far wake 
equations may be used. Schlichting [10] developed an analytical solution for 
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the velocity profile in the far wake region using a mixing length concept. In 
the far wake region the velocity profile is of the form 


V. - V 


10 [s 1-1/2 f 3/2" 2 

18 S C d tj 1 " b | 


( 8 ) 


where 8 is a constant that relates the mixing length to the wake width 
(8 - 0.18) , 


S is the streamwise distance from the vane trailing edge, 
y is the distance from the wake center, 
b is the wake half width: 


b - / 10 S(S C d g,) 1/2 an d 
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C d is a drag coefficient. 

Toe term C d £ is found from momentum consideration: 


C d * 


2 612 
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where 6^2 is the momentum thickness at the vane trailing edge and Vj/V* is the 
local velocity at the same location. The momentum thickness, 6j2» Is 
calculated by applying Truckenbrodt ' s energy-integral method [11] to calculate 
the turbulent boundary layer with a pressure gradient. 

After the wake profiles are calculated they are repositioned at the 
trailing edge of the casing vanes. The mass flow at this location is 
corrected for the wake blockage by integrating Equation (8) to give 
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OF POOR QUALITY 

The average velocity deficit over the wake is given as 

Q 'v w - v' 
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zu v - .max 

since for y - o at the wake centerline one has 
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V«» . max IS S Cj % 

The vane spacing is larger than the wake width so all velocities most be 
multiplied by the 
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to account for wake blockage, where L is the spacing for a particular vane. 

The effect of the calculated wakes on the potential flow field is shown 
in Figures 26 through 30 for the inlet casing vanes and the stator vanes. The 
wake half-width for each numbered vane calculated at the inducer inlet 
location is shown in Figure 31. Figure 32 shows the corresponding vane 
numbers. 

By making use of the streamf unction in the circumferential direction, the 
flow is calculated downstream. Thus the wakes from the inlet casing vanes are 
transferred to the smaller radius; and, because of the residual tangential 
velocity, are shifted to larger values of 0. In several cases the two waV.es 
combine to complicate the flow. For each of the wake distributions the 
massflow is integrated and checked. The curves will not all reach 180 degrees 
due to blockage of the boundary (fifth stator vane) at that point. 

In a similar manner the wakes are transferred to the inducer station in 


Figure 33. The circumferential shift is more obvious when comparing this 
figure to the previous ones due to the reduction in radius and streamwise 
distances traveled by the flow. 
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3. Experimental Verification of Analysis 

3,1 Rocketdyne Data 

The prediction of the Inlet flow to the Inducer Is a critical part of the 
damage analysis. The computational procedure has several key assumptions 
which need to be verified. The flow solution given In Figure 33 shows that 
the inducer inlet flow is non-uniform, producing large incidence variations 
which result In unsteady cavitation. This type of cavitation causes severe 
cavitation damage and high unsteady forces. 

Rocketdyne [12] conducted a series of experiments on the SSME HPOTP Inlet 
in air. For these experiments, a 5-hole probe [13] fabricated at ARL/PSU was 
utilized to obtain flow data in the inducer inlet plane. Total and static 
pressure, axial velocity and swirl angle were obtained as a function of 
circumferential angle at four radii. The axial velocity data are shown in 
Figure 34 and the swirl angle data are shown in Figure 35. 

A brief comparison between the predicted flow field given in Figure 33 
and the measured flow field given in Figure 34 shows similar trends. However 
the magnitude of the wakes are significantly different at the 2.23-lnch 
radius. It is proposed that these differences are due to separation occurring 
along the wall in this wake region. The trends in the axial velocity with the 
radius and the circumferential location are In good agreement. 

A better comparison is given in Figure 36 for the axial velocity and In 
Figure 37 for the swirl angle at the 1.977-inch radius. For this comparison 
the Rocketdyne data were Integrated to obtain an average velocity which then 
was used to normalize the data given In Figure 34. These results are in 


excellent agreement 
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3.2 Student Modal Data 

Three students from the Department of Aerospace Engineering designed, 
built and tested a 2-D model of the inlet casing. This work was part of an 
Undergraduate Thesis Program at The Pennsylvania State University. 

A full-scale model of the casing was constructed as shown in Figures 38 
and 39 with some modifications. The model was designed to be 2-D with uniform 
depth in the axial direction. All of the vanes were supported by pins so that 
their angle-of-attack could be varied. The shape of the throat area was the 
same as the inducer area on the turbopump. The inner throat and diffuser cone 
were made to rotate as a unit. This enabled a probe to be rotated 
circumferential./ at the Inducer inlet plane. A honeycomb was placed in the 
inlet to reduce incoming air turbulence. 

Results of using fluorescent mini-tufts on the casing walls and vanes for 
flow visualization showed that separation was occurring. To eliminate this 
flow separation on the casing, suction was applied at various points along the 
wall to keep the boundary layer attached. Flow separation was also observed 
on several vanes. 

Inducer inlet velocity data were obtained using a 5-hole probe. During 
these tests, the inlet flow velocity was approximately 175 ft/sec. The amount 
of suction required to maintain an attached casing boundary layer was set at 
a minimal value which was approximately 1% of the total flow. 

A comparison between data obtained by Rocketdyne with the 3-0 inlet 
casing and the student 2-D casing is given in Figure 40. Although the 
circumferential axial velocity profiles are in good agreement, the position 
and magnitude of the wakes are not. 
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4. Cavitation Analysis of Existing Design 

4.1 Observed Damage 

The locations of damage to the SSME-HPOTP are shown in Figures 3, 4, 5, 

6, 41, 42 and 43. Since the damage that currently limits the life of the pump 
is the deep pitting on the inducer casing, (Figures 4 and 5) the causes and 
recommendations for reducing this damage are of orimary consideration. The 
bands of cavitation damage on the casing (Figures 6 and 41), the light damage 
on the inducer blades (Figure 43), and the erosion of the inducer blade tips 
(Figure 3) will also be included in this analysis. The damage at the impeller 
exit (Figure 42) and possible causes for premature bearing failure will also 
be discussed. 

The deep pitting on the inducer casing as seen in Figures 4 and 5 is 
concentrated downstream of the stator vanes mainly on the counterswirl side 
and to a less extent on the preswirl side. The damage could not be due to 
cavitation on the statOV 'anes, since the local cavitation number is too high 
for cavitation to occur. For the flow over the stator vanes the cavitation 
number can be defined as: 



where P*, and V OT are the local values of static pressure and velocity. At the 
trailing edge of the stator vanes for the 109% power level* a = 4.7. This is 
much to high for any significant amount of cavitation. The damage must be due 
to other causes. 


*Note: 


All damage calculations will be for the 109% power level unless stated 
otherwise. 
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The local cavitation number on the inducer blade can be defined as 



where W*, is the local relative velocity of the blade section. For uniform 
flow and no preswirl at the inducer blade, a value of a = 0.108 was 
calculated. With a blade angle of 13.12° (r * 2.32 inches), the relative flow 
angle is 4.3° and a cavity length of approximately 33% of the chordlength 
should exist. The cavity length was approximated by the analysis discussed in 
Appendix A. At the hub a value of a = 0.318 was found and the cavity length 
was calculated to be only 7% of the chord length. The calculated cavitation 
pattern for this 'steady state' condition is superimposed upon the measured 
damage on the inducer blades and is shown in Figure 44. The trailing edge of 
the cavity at the tip corresponds to the region of damage on the blade. For 
this uniform flow condition the damage should be produced equally around the 
circumference of the casing. But how then could the damage be concentrated at 
certain areas on the casing? 

To answer this question the area where the damage is located provides 
part of the answer. In general the deep pits are downstream of the stator 
vanes. It is probable that, given the close proximity to the inducer leading 
edge, the wakes from the stators could be disturbing the flow. A decrease in 
axial velocity would increase the relative flow angle and thus the amount of 
cavitation. As the blades pass through these wakes some of the cavity is shed 
and the resulting collapses are creating the damage. For this reason the wake 
analysis was undertaken. 
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4.2 Cavitation Pattern Prediction 

As a first approximation the cavity length calculation will be performed 
assuming steady-state conditions locally in the wake. The local flow velocity 
at the inducer blade location can be calculated using the data given in Figure 
33. The axial and tangential velocity components and the flow angle were 
calculated through the wake for both preswirl and counterswirl sides of the 
pump. The cavitation number based on the local velocity on the blade was 
calculated using Equation (16). The cavity length was then calculated at 
points through the wakes using the analysis discussed in Appendix A. 

The cavity length calculations for the existing turbopump design are 
shown in Figures 44 and 45 and in Table 3. It would be expected that the most 
erosion should occur for both the preswirl and counterswirl sides between 24° 
and 51°, downstream of casing vane 3 and stator vane 15 (see Figure 32). As 
the inducer blade rotates circumferentially there would be a negligible cavity 
prior to entering the wake region. A large cavity would grow on the blade 
surface as the wake center is reached. Upon leaving the wake the cavity would 
be shed from the blade. Table 3 and Figures 44 and 45 also indicate a growth 
in the cavity as it passes through the wakes downstream of casing vanes 5 and 
6 and stator vane 16. For the preswirl condition the cavities exiting the 
wakes are longer than the cavities upon entering. The opposite is true for 
counterswirl. Thus, the volume of cavity shed should be larger on the 
counterswirl side of the pump. This correlates well with the observed damage 
patterns. 

4.3 Damage Rate Model 

The damage rate (D) is a function of a distribution of bubble collapse 
energies, the material response to the bubble collapse, and the duration of 
exposure. A general expression for the damage rate would take the form : 
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S 

D - P (Vi, a, o/o D , Ao(t) , To., Dn, h) • G(M, HT, D G ) • H(t/t 0 ) (17) j 

where for hydrodynamic considerations, 

) 

V» * reference velocity 

a * total gas content of the liquid F 

o/op * normalized cavitation number 

r 

| 

T,,,, * reference temperature 

D m = characteristic size 

i 

h = flow shape parameter i 

Ao(t) - time dependent cavitation number j 

for material response, 

M - test material (includes erosion resistance and corrosion 

resistance) 

I 

HT = heat treatments of the material ; 

Dq =» grain size/microstructure of the material (due to method of f 

fabrication, etc.) 
and for the exposure time 

t/t 0 » dimensionless exposure time. 

t 

i 

The total damage for cavitation damage sites would be . 

i 

n , i 

n T = Y Ditt (18) 

i = l 

and can be expressed as either a total volume loss or as a mean depth of 
penetration where the volume loss is divided by the damage area. 


* *js «... 
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Equatlon (17) for a particular case is very difficult to determine due to 
the many variables; however, for the HPOTP some assumptions can be made. 

First of all, one is only interested in comparing damage rates due to changing 
flow conditions. Thus, the material response enters as a constant and 
Equation (17) reduces to 

D = F(v£, o/a D , Aa(t) ) H(t/t Q ) . (19) 

For different flow conditions, the geometry of the pump also remains the same; 
however, the pump rpm and inflow conditions vary. 

The heavy cavitation damage that occurs on the housing has been shown to 
be caused by the change in cavity length as the inducer passes through wakes. 
Therefore, the damage rate is going to be different for each wake encountered 
by the inducer and Equation (19) reduces to a solution of Ao(t) which is 
approximately the change in cavitation volume. 

The change in cavity length as the inducer blade passes through a wake 
has been considered using a steady state analysis. It has been assumed that 
the cavity reaches a length corresponding to the incident flow at any point in 
the wake. If the wakes were large with respect to the blade chord this would 
be true. For the problem in question the wakes are narrow with respect to the 
chord and the time spent in a 20°-wide wake Is approximately 0.1 millisecond. 
The cavity would thus not have time to reach a steady state condition. The 
actual volume of the cavity shed can only be estimated for this condition. 

It will be assumed that (1) a steady state cavity is reached before the 
blade enters the wake and that a steady state cavity will exist some time 
after the blade passes through the wake, (2) a disturbance generated along 
the cavity wall will propagate with a velocity of W c = W^, ^ 1 + o which is the 
velocity along the free surface bounding a cavity, (3) the disturbance in the 
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cavity will follow a path corresponding to the flow at the leading edge where 
the disturbance was produced, and (4) the cavitation number is nearly constant 
through the wake. 

As a blade with a steady-state cavity starts to move through a wake, the 
blade will see a fluctuating flow angle due to the reduced axial and 
tangential velocity at each position. The change in incident flow angle will 
therefore create a disturbance at the leading edge of the cavity. The 
disturbance will propagate at the fluid velocity along the cavity wall. The 
magnitude and direction of the velocity vector at the leading edge is a 
function of the position in the wake. If the width of the wake is less than 
the length of the cavity on the blade as it enters the wake, then the blade 
leading edge will reach the other side of the wake before the initial 
disturbance to the cavity shape reaches the cavity trailing edge. This is 
important since for the wake cavity model there is no closure condition 
specified at the trailing edge. A closure condition could be specified but is 
not required for the cases that will be investigated on the inducer blade. 

The ’perturbed' cavity shape and shed volume can be calculated by the 
following procedure: 

(1) The inducer blade leading edge is placed at the exit point of the wake. 

(2) The flow angle and velocity at any point through the wake are used to 
calculate the cavitv contour at the sane point on the blade. 

(3) The volume of the shed cavity is 

AVol = Vol perturb - Vo ^exit * 
r ss 

Qualitatively, the results of such an analysis are shown in Figure 46. For 
case 2 in Figure 46, the steady-state cavity after passing through the wake is 
smaller than the cavity upon entering (typical counterswirl case for the wakes 
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from vanes 5 and 6). The shed cavity conies from two sources: the 

perturbation to the original cavity due to the wake, and the fact that the 
cavity volume must decrease even from the incident steady state condition. 

This occurs on the counterswirl side since the axial velocity is increasing as 
the inducer blades move from 180 8 to 0° , For case 3 in Figure 46 (wakes from 
vanes 5 and 6 on the preswirl side) the opposite is true. The cavity leaving 
the wake is typically larger than the cavity on the blade when entering the 
wake. Some of the volume in the 'perturbed' cavity must be used to 'fill' the 
cavity leaving the wake. If the wake is not very 'severe', there may be no 
shed cavity at all. As can be seen from Figures 6 and 41 for the casing 
damage, the erosion on the counterswirl side is much worse than on the 
preswirl side for the positions behind vanes 5 and 6 and stator vane 15. 

For both preswirl and counterswirl sides the heaviest pitting was 
downstream of casing vane 3 and stator vane 15. The steady state cavities 
entering and leaving the wake are small (less than 0.1 l/c) for both preswirl 
and counterswirl. Thus nearly all the cavity volume generated while the blade 
was in the wake would be shed. The analysis for this wake can be applied if 
the incident and exit cavities are neglected (this is done because no closure 
condition is used for the cavity wake model) . This cavitation should be very 
severe since the entire 'perturbed' cavity is shed. The amount of damage can 
be approximated by the volume of the cavity shed. 

Emphasis has been on damage from unsteady cavitation created when the 
inducer enters a wake. An analysis is presented in Appendix B which clearly 
shows that this type of damage is at least an order of magnitude greater than 
its steady capitation counterpart. 
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4.4 Cavitation Damage at the Impeller Exit and Premature Bearing Failure 

The cavitation damage at the impeller exit Is due to the difference in 
pressure across the blade at the trailing edge. The high pressure fluid on 
the pressure side will flow towards the low pressure region on the suction 
side. Because of the sharp trailing edge geometry this flow will separate. 

Even though the exit pressure is almost 5000 psia, the local cavitation 
number at this location is a * 0.818 since the impeller tip speed is 878 
ft/sec. This means that cavitation would occur in a separated region. 

Breaking the sharp edge at the impeller blade trailing edge by rounding the 
pressure surface should reduce this problem with a minimum loss in 
performance. 

The premature bearing failure is likely due to the unsteady forces 
generated by the unsteady cavitation on the inducer blades. The proposed 
design change should also reduce this problem. 

5 . Recommendations for Cavitation Damage Reduction 

5.1 Casing Vane Modifications 

The analysis of the existing stationary vanes in the HPOTP clearly 
indicates that the vanes must be altered to improve the turbopump's cavitation 
performance and hence reduce cavitation damage. Modifications will strive to 
minimize the fluid turning at the inducer, equalize the massflow around the 
circumference of the inlet, and reduce the sensitivity to flow incidence. The 
last of these criteria will be considered first. 

There is a bonafied reason to use airfoil shapes in the redesign of the 
inlet casing vanes. Regardless of the camber or amount of turning that a vane 
is to accomplish, typical vane shapes always employ a rounded leading edge and 
a sharp trailing edge. Due to the uncertainty in the exact angle at which the 
flow will encounter the vane at the leading edge, a rounded surface will allow 
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the flow to oass beyond the leading edge without flow separation on either the 
pressure or suction surface. The trailing edge Is made as thin as possible to 
minimize the possibility of flow separation. 

A otaiiJard NACA 65—0012 series vane shape was used as a base design for 
the new stationary vanes. The camber offset was increased for each vane until 
the desired amount of vane turning was produced. For this analysis the 
trailing edge thickness was set equal to zero. The thickness was 
added normal to either side of the camber line up to the maximum thickness at 
30 percent of the chord. At this point it was then added to the suction side 
in an exponential manner until at the trailing edge all of the thickness was 
added to the suction surface. 

The Douglas Neumann Cascade Program solves for the flow at any location 
in the cascade. Although the flow at the trailing edge must be in that 
direction, the inviscid flow between vanes can still exhibit a substantial 
degree of turning different from that indicated by the adjacent vanes. Thus 
some overturning may be required to direct the flow radially at the exit of 
the stationary vanes on the average. 

Several redesigns were considered. Although the analytical results show 
that a substantial improvement in the flow field can be obtained, the vanes 
should not be considered a final design. As will be pointed out, the casing 
boundary plays a significant role in the guidance of the flow, and a new 
design should incorporate a redesign of the casing as well. 

The stator vanes serve no useful purpose hydrodynaraically . While they 
might be essential structurally, their presence in the flow field only serves 
to aggrevate the situation at the inducer because of their close proximity. 

The current redesign employs six vanes instead of four but makes them 
one-third shorter. The actual number and size of these vanes should be 
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dictated by structural considerations. The bottom stator vane should be 
incorporated into the redesign of the inlet casing boundary. While the new 
stator vanes will locally divert the flow, they will have no impact on the 
flow at the inducer and could readily be removed from the design 
considerations. 

The inlet casing vanes were increased in number from 13 to 15. This will 
serve to maintain the force level on each vane since more turning will now be 
expected of the vane row to eliminate the tangential component of the 
velocity. They were spaced equally at their trailing edge (radius of 4 
inches) between the top and the casing boundary at the bottom. The amount of 
turning was calculated to allow the flow to enter the vane row with no 
appreciable Incidence and to exit radially. In some cases fluid turning 
angles close to 90 degrees was indicated. This is typical for inlets of this 
type where the flow is first directed tangentially by the casing near the 
upper inlet and then must be turned radially, effectively increasing the fluid 
turning before reducing it. Again, this emphasizes the important roll that 
the design of the casing boundary plays. 

The vanes were uniformly decreased in chord length from the top casing 
vane, number 2 to the bottom vane, number 9. Vertical symmetry was 
employed to balance the side forces on the shaft. Obviously the downward 
force will still be substantial. The inviscid flow was calculated to produce 
radial flow at the radii corresponding to the inducer station, thereby 
eliminating any tangential velocity. 

A typical vane shape of 45 points is illustrated in Figure 47. This vane 
is nondimensionalized by the chord and the dotted camber line illustrates how 
the thickness near the trailing edge is shifted towards the suction surface. 

In all, 22 bodies (vanes) are now used in the cascade with a total of 1092 
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points. The redesigned cascade took 2 1/2 hours to run on the computer, fhe 
flew solutions and new vane shapes are illustrated in Figure 48. The stators 
are now cut at the trailing edge at the same radius from inner to outer 
streamline and the wake flow produced by these two vane rows is shown in 
Figure 49. These wakes will be more diffused when they reach the inducer due 
to the increased distance between the stator trailing edge and the inducer 
leading edge. 

The flow field at the inducer leading edge is shown in Figure 50. With 
the reduction of the tangential velocity the flow becomes more uniform at the 
inducer and the circumferential massflow distribution is nearly constant. 

Note that the casing boundary is still responsible for some fluid turning 
towards the bottom of the inlet at large values of 0. This turning can only 
be removed by redesigning the bottom of the inlet casing. 

As can be noted in Figure 50, the deepest wakes are produced by the six 
stator vanes. It appears that due to the location of these vanes (in close 
proximity to the inducer leading edge) the wakes will always be the worst and 
perhaps be the controlling factor for cavitation damage. It is recommended 
that these vanes be removed if possible. 

The va u sh pes in Figure 48 were designed to produce radial flow by an 
inviscid flow analysis. Due to the viscous boundary layer, additional turning 
will have to be placed in the vanes to achieve the same turning as predicted 
by the inviscid solution. This is accomplished by adding the flow deviation 
angle to the blade camber angle. The leading and trailing edges were held in 
place (thus the chord remains constant) and the camber increased to include 
the fluid deviation angle. The camber offset distribution was altered to 
place most of the additional camber towards the trailng edge so as not to 
disturb the flow and hence the incidence angle at the vane leading edge. The 
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newly designed vat ^s with offset a rn illustrated In Figure 51. A typical vane 
Is Illustrated In Figure 51 and can be compared with Figure 47 to see the 
effects on the camber of the fluid deviation. Coordinate data for the 
redesigned vanes are given In Appendix C. 

If It is required to thicken the trailing edges of the redesigned vanes, 
this thickness can easily be added to the suction surface as indicated by the 
dotted line in Figure 47. In this way the flow is still guided by the 
pressure surface and the vanes are thick enough to be cast. 

The vane offsets were estimates from a deviation analysis similar to that 
discussed in Section 2.2. Using cascade data for the vanes in the casing Is 
an approximation which has a high risk. The proposed vanes do more flow 
turning than the existing vanes and hence larger offsets are required. Due to 
the lack of available data for vanes in the accelerating flow created by the 
casing, it is strongly recommended that air tests be conducted to verify or 
adjust the vane geometry to meet the rquired flow performance. 

5.2 Cavitation Damage Analysis 

The cavity length calculations, based on the worst wakes, for the new 
designs with and without stator vanes are shown In Figure 53 and in Table 4, 
There are a number of improvements for both redesigns over the existing 
turbopump. These can be summarized as follows: 

(1) The flow velocity Is more uniform, thus the cavity lengths entering and 
leaving the wake are nearly equal. 

(2) The velocity deficit is reduced significantly for the redesign with no 
stator vanes. 

(3) The widths of the wakes are reduced. 
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The expected reduction in the rate of cavitation damage la proportional 
to the change In shed cavity volume. For the worst wakes, comparing the 
existing design to the redesign without stators. 


^old (shed cavity volume)^ 


(shed cavity volume) 


- 21.1 


( 20 ) 


'new 


new 


6. Summary 

A model has been developed that correlates the severe cavitation damage 
on the high pressure oxidizer turbopump (HFOTP) casing with the inducer 
unsteady cavitation. A computational procedure was used to analyze the flow 
through the inlet casing including the prediction of wakes downstream of the 
casing vanes. The predicted flow at the Inducer Inlet compares very favorably 
with the experimental data obtained by Rocketdvne. 

The analysis of the HPOTP clearly shows that the inlet casing vanes must 
be altered to significantly reduce the cavitation damage rate. However, any 
redesign should reduce the vane sensitivity to flow incidence, equalize the 
distribution of the mass flow around the circumference of the inlet, minimize 
the swirl at the inducer Inlet and, if possible, eliminate the stator vanes In 
front of the Inducer. 

Using the damage analysis presented, the damage rate for the proposed 
redesign is reduced by a factor of 21. The proposed redesign uses a 
conventional vane shape and the vanes are cambered more than the original 
vanes to significantly reduce the swirl component of the velocity. In 
addition, the Inlet casing vanes were Increased In number from 13 to 15 to 
reduce the forces on the vanes at the bottom of the casing. 


The following conclusions can be made: 
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(1) The cavitation damage problem on the HPOTP results from the unsteady 
Interaction of the inlet casing and stator vane wakes and t|ie cavity flow 
over the Inducer blades. 

(2) The existing vanes produce separated flow and large wakes which combine 
to give severe velocity deficits at the leading edge of the Inducer. 

Also, a significant amount of fluid turning exists downstream of the 
vanes. 

(1) The stator vanes should be removed, if possible, due to their close 
proximity to the inducer blades. 

(4) A redesign of the inlet casing vanes is required to greatly improve the 
cavitation damage rate occurring on the housing. 

(5) A complete redesign should address the shape of the casing boundary as 
well. 

(6) The redesign should also reduce unsteady cavitation and bearing forces. 

An analysis was conducted which shows that the non-uniform flow entering 
the inducer causes the severe cavitation damage that exists in the HPOTP,. The 
redesign of the stationary vanes is an initial attempt to significantly reduce 
the cavitation damage rate. It is strongly recommendtd that a series of tests 
be performed on any redesigned inlet casing to verify that the improvement in 
the cavitation damage rate will be met. These tests should be conducted along 
with continued efforts at refining the computational analysis. 
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Table 1 

SPACE SHUTTLE MAIN ENGINE CHARACTERISTICS 


Engine Operation 
Thrust 

Vacuum 

Chamber Pressure 
Area Ratio 

Specific Impulse (Nom) 
Vacuum 

Mixture Ratio 
Life 

Specification Dry Weight 


MPL - RPL - FPL 

360K - 470K. - 512K 

1957 PSIA - 3012 PSIA - 3283 PSIA 

77.5 

455.2 

6.0 

7.5 Hrs 
55 Starts 

6957 Lbs 
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Table 2 

HPOTP OPERATING CONDITIONS 


Power Level 

MPL 

RPL 

FPL 

Percent Thrust Level 

65 

100 

109 

Pump Inlet Pressure, PSIA 

3 4 

377 

390 

Pump Inlet Temperature, °R 

*68 

170 

171 

Pump Inlet Density, lb/ft^ 

70.5 

70.2 

70.1 

Pump Inlet Flowrate, lb/sec 

712 

1068 

1160 

Pump Dlsch Pressure, PSIA 

2513 

4290 

4790 

Pump Disch Temperature, °R 

180 

192 

196 

Pump Disch Density, lb/ft^ 

70.6 

70.3 

70.2 

Pump Head Rise, ft 

4474 

7915 

8893 

Pump Speed, RPM 

20298 

28349 

30367 

Pump Power, HP 

8385 

22571 

27663 
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Table 3 

Calculation of Cavity Lengths for the Existing Design 

(Worst Wakes) 

R = 2.37 Inches 


Position Flow Angle Cavity Length 

9, Degrees o, Degrees Chord Length 


Preswirl Side 

Entering Wake 3 

24.3 

-1.1 

0.08 

Pressure Side 

In Wake 3 

35.9 

6.4 

0.40 

Suction 

Side 

Leaving Wake 3 

50.7 

0.2 

0.02 

Suction 

Side 

Entering Wake 5 

81.9 

2.2 

0.20 

Suction 

Side 

In Wake 5 

88.7 

7.0 

0.43 

Suction 

Side 

Leaving Wake 5 

95.4 

3.2 

0.27 

Suction 

Side 

Entering Wake 6 

95.4 

3.2 

0.27 

Suction 

Side 

In Wake 6 

110.4 

8.1 

0.46 

Suction 

Side 

Leaving Wake 6 

121.9 

4.1 

0.34 

Suction 

Side 

Counterswirl Side 

Entering Wake 3 

50.7 

1.3 

0.08 

Suction 

Side 

In Wake 3 

35.9 

6.5 

0.40 

Suction 

Side 

Leaving Wake 3 

24.3 

0.1 

0.00 

Suction 

Side 

Entering Wake 5 

95.4 

4.3 

0.35 

Suction 

Side 

In Wake 5 

88.7 

7.5 

0.44 

Suction 

Side 

Leaving Wake 5 

81.9 

3.6 

0.30 

Suction 

Side 

Entering Wake 6 

121.9 

4.7 

0.36 

Suction 

Side 

In Wake 6 

110.4 

8.5 

0.47 

Suction 

Side 

Leaving Wake 6 

95.4 

4.3 

0.35 

Suction 

Side 
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Table 4 

Calculation of Cavity Lengths for the New Designs 
(Worst Wakes) 

R => 2.37 Inches 


Without 

Stators 


With 

Stators 


Position 
9, Degrees 

Flow Angle 
a. Degrees 

Cavity Length 
Chord Length 

51.2 

3.6 

0.30 

55.6 

6.0 

0.39 

59.6 

3.6 

0.30 

80.8 

3. 8 

0.32 

85 

7.0 

0.43 

87.5 

3.8 

0.32 


POWERHEAD COMPONENT ARRANGEMENT 



OXIDIZER 

PREBURNER 




HIGH PRESSURE OXYGEN TURBOPUMP 
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Fitrnr.; 2. Schematic of the HPOTP (RocketHyne) 
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LIGHT 
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APPENDIX A 

Steady-State Cavity Calculation Procedure 

A method developed by Stripling and Acosta [14] was used to calculate the 
size and shape of the cavity on inducer blades. The flow was modeled for a 
helical inducer by a free-streamline calculation of the flow through a 
semi-infinite cascade. The flow was assumed to be two-dimensional, 
irrotational and inviscid. A wake or dissipation model simulates the blockage 
that occurs due to the mixing at the cavity trailing edge which is the 
condition emphasized when the cavity is less than the chord length. The flow 
is then confined within the cascade blade passage and the chord assumed to be 
infinitely long. 

The free-streamline flow is mapped in the hodograph plane to simplify the 
cavity calculations. A schematic of the flow in the physical and hodograph 
planes is presented in Figure A1 . The hodograph variable is defined as; 

C = 4— = u - iv (A.l) 

s dz 

where u and v are the x and y velocity components in the physical plane and 
z = x + iy. 


+ i (p 


(A. 2) 


is the complex potential of the flow with <j> being the velocity potential and 
the stream function. The physical coordinates of the flow are found by 
determining the functional relationship between F and 5 and integrating A.l; 


t dF 

z ■ J — + const 


(A. 3) 
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To solve A. 3 the following boundary conditions in the hodograph plane must 
hold: 

(1) The real axis and |c| * W c are streamlines 

(2) ; * 0 is a stagnation point, therefore 

f|( C -0) -0 

(3) Streamlines in the hodograph plane, near the origin corresponding to the 
flow at infinity, are uniform and inclined at an angle a to the axis. 

From Bernoulli's equation the velocity along the cavity wall W c is: 

W c - tf* j 1 + a (A. 4) 

where is the upstream velocity and a is the cavitation number given as: 



The coordinates of the flow are then found from A. 3. The coordinates of the 
cavity length and height are found from: 

z ■ x + iy 
3 c + ih 

“ z(£ 3 W c ) - z( j; = - W c ) 
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APPENDIX B 

Unsteady Cavitation and Erosion 

The destructive action caused by cavitation has been a practical problem for 
many years. Erosion of ship propellers and hydraulic turbine blades can cause 
a loss in performance leading to eventual costly maintenance or replacement of 
the damaged parts. The damage is caused by the collapse of the vaporous 
cavities as they move into a high pressure region. Upon collapse, enormous 
pressures are generated. If the cavity collapse occurs near a boundary, 
damage can occur. 

Accurate prediction of cavitaiton damage is very difficult. After nearly 
100 years of research in cavitation, the problem of calculating the amount of 
erosion for a given prototype Is as yet unsolved. The lack of success is 
because cavitation damage Involves both the fluid and solid mechanics and thus 
is Inherently difficult. 

The attack caused by the collapse of one cavitation bubble occurs over a 
very small area (or the order of hundredths of a square millimeter) and in a 
very short time interval (measured in microseconds). Both the hydrodynamic 
aspects of the cavity flow and the material response to the impact loading 
caused by the cavity collapse trust be considered. The problem is further 
complicated by possible interactions between cavitation and corrosion. 

In analyzing the Inducer blade and casing damage it was assumed that a 
steady sheet of cavitation would exist on the blade in uniform flow. Any 
cavity length fluctuations arise from the blade passing through the wakes. 

This sheet type of cavitation is often observed near the blade tip on surface 
ship propellers. Most conventional surface ship propeller blades have low 
thickness to chord length ratios and low camber. The typical blade sections 
at the tip are similar to the inducer blade profile for the oxidizer pump. 


It has been reported that when a surface ship propeller blade passes 
through the wake caused by the hull, severe noise, vibration and structural 
damage (blades and hull) can result. In a recent paper presented at the Joint 
Symposium on Design and Operation of Fluid Machinery, June 1978, Cumming, 
et.al. [15] measured the effect of wakes on the cavitation damage rate for 
some surface ship propellers. The results for two tests, one with and one 
without an upstream wake, are presented in Table Bl. The mean velocity, test 
duration and cavitation number were approximately equal for both tests. Since 
the flows were similar except for the wake any change In damage rate should be 
due to the unsteady flow. The total pitting rate was approximately 2.48 times 
higher for the flow with the wake. The wake was approximately 60° wide. If 
the difference in the pitting rate were due solely to the wake then the 
pitting rate in the wake is nearly 8.9 times higher. While the pitting rate 
is one measure of the damage the size of the individual pits should also be 
considered. 

For an approximation the pits are assumed to be conical sections with a 
constant diameter to depth ratio. The damage data is reduced to a total 
volume of erosion per unit area per hour. The volume erosion rate is 140 
times higher for the wake region than in uniform flow. Thus, decreasing the 
velocity deficit and wake width can greatly reduce the damage rate. The 
assumption that the pits have a constant diameter to depth ratio probably is 
conservative. It was shown by Stinebring [4] for a cavitation damage study 
in the incubation zone that generally, the depth/diameter ratio increases with 
increasing diameter. 

Another study of surface ship propellers by Wilson et.al. [16] for Naval 
Auxiliary Oilers (AO-177 class) studied the vibration and damage problem due 
to wakes. The ships had severe vibration transmited through the hull and 
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damage on the propeller blade trailing edges near the tip* They attributed 
the cause of the problem to the unsteady cavitation on the blades as they 
passed through the hull wake* The cavity would grow while in the wake, then 
break off and collapse against the trailing edge when leaving the wake. 

Adding a special fin upstream of the propeller reduced the wake deficit which 
reduced vibration and damage. Observations also showed that In addition to 
reducing the volume of cavitation, the cavity also collapsed in the free 
stream away from the blade. 

These are just two examples of the increased damage due to the unsteady 
cavitation produced on blsdec passing through wakes. 
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Table Bl 


Run U5 a ■ 3.44 

Run W1 o ■ 3.51 


Pitting Rate ■ 540 pits/cm 2 /hr 
Pitting Rate - 1340 pits/cmfyhr 


Size Distribution of Pits 
(X of Pits Within Size Range) 


2R 


Diameter 

(ram) 

Run U5 

Run 1 

< 3.0 

36 

9 

3. 0-4.1 

47 

19 

4. 2-5. 3 

11 

12 

5. 4-6. 4 

3 

13 

6. 5-7. 4 

3 

12 

7. 5-8. 4 


12 

8. 5-9. 5 


10 

9.6-10.6 


3 

10.7-11.7 


3 

11.8-12.9 


2 

13.0-14.0 


l 

14.1-15.0 


l 

15.1-16.1 


2 

16.2-17.3 


1 


Average Over 
One Revolution 


Average Over 
Wake 


V - 1/3 R 2 (C(R)) 

- (C/3) ir R 3 


Volume Los 8 
Rate 


(Wl) 89.2c 
(U5) " 3.7c 


- 24.3 


Volume Loss 
Rate 


Wl - U5 360* 


W~ 


Pitting 

Rate 


Wl - U5 // 360 #N 
U5 60° 


- - — * 
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APPENDIX C 

Coordinate Data for the Redesigned Inlet Casing Vanes 

The following sets of coordinates can be used to construct the newly 
designed vanes. All dimensions are in inches and the vane shapes are numbered 
as in Figure 48. The x-axis is positive to the right of the page and the 
y-axis is positive to the top of the page. 
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